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Abstract 

A theory for surface transitions in the presence of a disordered pinning poten- 
tial is presented. Arbitrary disorder correlations are treated in the framework 
of a dynamical functional renormalization group. The roughening transition, 
where surface roughness and mobility behave discontinuously, is shown to turn 
smoothly into the continuous superroughening transition, when the range of 
disorder correlations is decreased. Implications for random-field XY-models 

and vortex glasses are discussed. 
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The shape of a fluctuating crystal surface is affected by the pinning potential provided by 
the crystal planes. The interface exhibits a roughening transition, if these planes are perfect, 
and it exhibits a superroughening transition, when these planes are strongly disordered. The 
roughening transition is well understood theoretically and verified in different experimental 
situations (see e.g. [|I],0 and references therein). Above the roughening temperature T R the 
interface is thermally rough, whereas below Tr the interface becomes smooth, since it locks 
into the pinning potential. However, for quenched bulk disorder |§ and substrate disorder 
@H it was shown, that the transition has a different nature. This so-called superroughening 
transition occurs at Tsr = Tr/2. Above Tsr the interface is thermally rough again, but 
below Tsr the disordered pinning potential increases the roughness of the interface. In this 
Letter we address the question: how are the roughening transition and the superroughening 
transition related? What happens "in-between" , if we switch from a pinning potential with 
long-range correlations to short-range correlations? 

As shown below, new physics emerges for logarithmic substrate roughness. This can be 
realized, if the crystal is grown on a logarithmically rough substrate. Such a substrate could 
be generated by quenching it from a temperature above its roughening transition. In order 
to have a stable substrate surface, the roughening temperature of the substrate should be 
higher than that of the crystal. 

Our analysis is also relevant for a AF-model in the presence of a random field, as long 
as vortices can be neglected ||. Different types of sources of the random field may have 
different multipolar character, which lead to a power-law decay of correlations. 

Finally, the consideration applies to two-dimensional vortex glasses [|]||. Their pinning 
can be caused by local suppressions of the condensate density, which also relaxes with power 
laws and leads in turn to a power-law decay of correlations. This power varies with the 
dimensionality of the defects. 

At present there is a controversial discussion about the actual roughness of the super- 
rough phase. Renormalization group (RG) calculations [|||, variational calculations || and 
simulations [ lOj give an inconsistent picture. In view of this controversy we focus on topics 



consistently described by these different approaches: the phase boundary and the nature of 
the phase transition. The present study is based on a functional dynamical renormalization 
group applying to the limit of very weak pinning. 

The interface profile is described by as a function of time t and the two-dimensional 
lateral coordinate r. Its kinetics is captured by the overdamped equation of motion 

m' l d t (p{t, r) = KV 2 <p{t, r) - x sin [<j>(t, r) - d(r)] + 

+F + £{t,r). (1) 

Here m is a mobility, K is the surface stiffness, % * s the amplitude of the periodic pinning 
potential, which has a quenched disorder phase d representing deformed crystal planes, F is 
a driving force per unit area, and £ is thermal noise at temperature T. We implicitly suppose 
a cutoff A for the wave- vector of shape fluctuations, which corresponds to an area 4n/A 2 per 
degree of freedom in the lateral plane. Since we assume regularization in momentum space, 
we consider <fr and d to be defined continuously on the whole lateral plane (see Fig. [I]) and 
not just only on a lattice of points r |[L1|| . 

The disorder field is characterized by zero mean and a difference correlation 



[d(r) — d(r')] 2 = 2A(r — r'), which for simplicity is taken to be isotropic. For the following 
it is convenient to introduce the correlation (g := \m 2 x 2 ) 



7o (r) := go e*W*M(o)] = g Q e" A « . (2) 

In terms of this function, a perfect crystal exhibiting the roughening transition has 7o(r) = g. 
In the studies of the superroughening transition 7o(r) was supposed to decay rapidly 

for large r. Eq. (fj) shows, that this case corresponds to a very rough substrate with a rapidly 
increasing function A(r). Here we allow for a general, possibly slowly decaying correlation 
7o(r). 

In constructing the RG, we choose an approach different from previous treatments. Since 
we do not want to specialize to a special correlation, we have to perform a functional RG. 
This can be achieved in the formalism of Martin, Siggia, and Rose |L2|| , which requires the 



introduction of an additional field (p. In this formulation the disorder average can easily be 
performed. The generating functional then reads 



V(j)V(j) exp[A ] (3) 
with the action Aq = Aq + Aq composed of free part and disorder part 

4 0) = Jdt d 2 r {j o (t,r)0 o (t,r) + ^ $(*,r) - 

-0 o (t,r) ^fafar) - K V 2 (j)o(t,r) }, 
Af ] = Jdt dt' d 2 r d 2 r' i 7o (r - r') x 

x0 o (t, r)0 o (t', r') cos[0 o (t, r) - o (t', r')]- (4) 

We introduced the unrenormalized quantities O — 4>, $0 — 2mT, /i = 1, n Q = mK, and 
Jo = mF. 

In the absence of the pinning potential (x = and 70 = 0) one has correlations 

C7<°>(f,r) :=^[0(t,r)-0(O,O)] 2 )W 

= nC i0) (/iKA 2 t,Ar), 

C (0) (r,p) = 2/ 1 ^[l-J (xp)e-^l] 
jo x L J 

ps lnmax(l, |r|, p 2 ). (5) 

Jo denotes the Bessel function and n := fi^/(8irK) = T/(AttK) is the roughness coefficient 
of the unpinned interface. 

We perform a momentum-shell RG [13] on model ([3]) integrating out all modes with 



wave vectors Ae~ dl < k < A of both fields <fi and <fi. This calculation is performed only to 
first order in 7 without restricting its functional form. In order to restore the original value 
of the cutoff, we then rescale r — > e dl r, and simultaneously t — > e 2dl t, <pi — > e~ 2dl <fti + di, and 
(f>i — > (pi+di according to the scale invariance in absence of pinning. 
This procedure easily gives the flow equation for the correlator 

$ T ,(r) = (4 + r.V-2n,)7j(r), (6) 



which is exact in this order. On the same level, there is no renormalization of /i, and k. 
However, new interaction terms are generated, which are not present in the unrenormalized 
action. In order to keep the theory simple, i.e. the number of parameters small, we replace 
these terms approximately by contributions to the flow of /i, and k. Following Nozieres and 
Gallet [0, we thereby take the coupling between the smooth modes (small k) into account 
in order to avoid artifacts related to the special choice of the cutoff procedure. The resulting 
flow equations are (p = Ar) 

1 



diKi 

dm 



87TfC;A 4 

2vmfA 4 

-fa* 



d 2 p p 2 ll (p/A)A(n u p) 
d 2 p ll (p/A)B(n h p) 



diJi = 2Ji 



(7) 



with the definitions 



ni :-- 
A(n,p) :-- 
B(n,p) :-- 



8ttki 
J (p) e 



-nC(°)(0,p) 



J (p) / dr e" 
Jo 



-nC(°)(r,p) 



(8) 



Under renormalization the effective stiffness of the interface is increased, diKi > 0, and 
the roughness amplitude is decreased, dpni = — (n// Ki)diKi < 0. As long as disorder is renor- 
malized to zero, it will have the effect to reduce the equal-time roughness given by n^, since 
C(t, r) ~ n oo C^ \p oo K oo A 2 t, Ar). The invariance di{p,i'&i) = expresses the preservation of 
the fluctuation-dissipation theorem under renormalization. This allows a stepwise solution 
of the flow equations: first for 7 and k, since their flow equations are closed. These quan- 
tities completely determine the static properties of the system. Afterwards one can solve 
for 1? or /i, which affect only the dynamics. Pinning always reduces mobility, dipi < 0, even 
if disorder is renormalized to zero on large length scales. The driving force leads to a drift 
velocity (0) ~ \x\*mF . In the driven case renormalization has to be stopped at a value I* of 



the flow parameter determined by mF = e~ 21 * K;*A 2 since a finite drift velocity washes out 
the effect of pinning. 

Now we discuss the phase boundary separating the regions, where disorder changes from 
being irrelevant to being relevant. The flow equation for 7 can directly be integrated: 



7z(r) = exp 



2 /' dl'{n v 
Jo 



7o(e'r). (9) 



In two extreme cases the relevance of disorder is obvious. The first case 70 (r) = go represents 
perfect crystal planes. From eq. (||) we find the boundary at = 2 which is the usual 
roughening transition. In the limit \ ~ * of a weak potential the transition temperature 
is Tp°' ) = 8ttK. For finite x the transition temperature Tr will be higher than \ since 
^00 < no- In the second case 70 (r) = go5(Ar) of extremely deformed crystal planes we 
find the boundary at = 1. The superroughening transition occurs at a temperature 
T = Tsr = 4nK independently of disorder strength, since here always rioo = no- This 
universality of the transition temperature is special for a infinitely sharp 70 (r). A finite 
width leads to d\ni < and hence an increased transition temperature Tg^ > 4nK. Such 



a shift has been observed recently in simulations of superroughening \Ti\. In a comparison 
with such data, the lateral lattice spacing should be taken as width. 

The decisive interpolation between theses cases of very broad and very sharp 7(1") is 
provided when 7 (p/A) w Cop~ 2a for large r with some amplitude Co > and some ex- 
ponent a > 0. This means, that the crystal planes themselves have logarithmic rough- 
ness. Renormalization preserves the asymptotic exponent a, but the amplitude is changed: 
q oc exp[— 2 Jo dl'{ni> + a — 2)}cq. Thus we expect pinning to become relevant for the large- 
scale behavior of the interface at rioo = 2 — a. 

To be more specific, we solve the flow equations (0), after eq. @ has been plugged in. 
For weak pinning, an eventual transition is related to large values of the flow parameter, 
where the flow equations simplify: 



dmi = -y\ 



2 
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yf 



A < e -2j l dV { n v+a -2) for a<2 

A> e - 2 fo dl ' n i' fora>2 



B < e -2^dl'(n ll+ a-2) for Q < j 

dainty « 4 , (10) 

B>e - 2 / ^(n ; ,-l) fora>1 

The various coefficients like A < follow from the original equations (|7|) using suitable approx- 
imations for large I, which are not displayed here since they are clumsy and straightforward 



15fl . These coefficients are positive and well-defined at least for n > 5/4 — a. This region 
includes the transitions under discussion. As long as we are not too close to that border, we 
may ignore the implicit dependence of the coefficients on n. Due to the ambivalent asymp- 
totic behavior of flow equations (|i~0|) one has to distinguish three regimes depending on a 
displayed in Fig. |[ 

Regime < a < 1. The flow trajectories in the (n, y)-plane have hyperbolic structure 
(ni + a — 2) 2 — yf =const. like in the Kosterlitz-Thouless transition ||16|| . In addition, 
$i = i9o ex P[ — {A K / B < )(rii — no)]- Disorder is irrelevant as long as n + a — 2 > y since 
then t/oo = and > 2 — a. For weak disorder strength, the boundary = 2 — a 
corresponds to a temperature Tr(«) ^ (2 — a)A7cK. Approaching the transition, and 
$00 remain finite, as well as the effective mobility /i^ = /io^o/^oo- Thus for < a < 1 
the transition is qualitatively similar to the roughening case a = 0. In the low-temperature 
range, the perturbative flow equations loose their validity, since the strength of pinning 
increases. However, as the effective amplitude q of 7/ diverges as discussed above, we expect 
the interface to lock into the pinning potential. This implies, that the equal-time roughness is 
C(0,r) ~ aln(A 2 r 2 ) below the transition, whereas it is C(0,r) ~ nooln(A 2 r 2 ) above. Hence 
the transition is discontinuous. Following a — > 1~ along the transition, the magnitude of 
the jump disappears, since for a = the logarithmic roughness of the crystal planes and 
interface coincide. In a similar way the (linear) mobility behaves discontinuously at the 
transition since it remains finite approaching the transition from above and it is zero for the 
locked-in interface. The jump of the roughness at the transition from = 2 — a to a is 
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universal, whereas the jump of mobility depends on all parameters. 

Proceeding to the regime 1 < a < 2, only the flow equation for *& changes. It exhibits now 
a continuous divergence, if approaches the value 1 from above, where Tsr(o;) ^ 4ttK. 
Thus the interface becomes immobile (/Xqq = 0) at that point, although it does not lock in, 
what happens only at the lower temperature T R (a) given before. Thus our flow equations 
suggest a sequence of two transitions: a dynamic one and a static one. However, due to 
their perturbative nature, we are not really in position to characterize the low-temperature 
phase n < 1. Possibly the lock-in transition for 1 < a < 2 might be an artifact of the first 
order calculation. If the transition exists, the interface is rougher below the transition than 
above. 

In the third regime 2 < a we also have to use the second flow equation of n. Nevertheless, 
the flow remains hyperbolic with nf — yf =const and every n > is renormalized only by 
a small amount. Thus disorder is always irrelevant for n, a lock-in transition becomes 
impossible. The dynamic transition is the same as in the regime 1 < a < 2. Again, we can 
not seriously say something for n < 1, except that the linear mobility should continue to be 
zero. The controversy mentioned in the introduction is precisely about this phase, since the 
commonly used 7(1") decaying faster than any power belong to a = 00. However, there is 
agreement that below the superroughening transition the interface is at least as rough as at 
the transition. This fact is not captured by the first order flow equations. 

To sum up, we have performed a dynamical RG for a disordered sine-Gordon model 
in the presence of weak pinning with arbitrary disorder correlations. Thereby roughening 
and superroughening transition could be merged for a correlation 70 (r) ~ r~ 2a . Disorder is 
irrelevant for the interface at temperatures T ^ max(2 — a, 1)4itK. In the range < a < 1 
the transition is of roughening type with a discontinuous jump of roughness and mobility. 
The jump of roughness is universal. This discontinuity decreases, as a = 1 is approached. For 
a > 1 the transition is of the continuous superroughening type. The transition temperature 
is non-universal in all regimes, as long as disorder correlations have a finite range. 

A similar interpolation between transitions has been found in a study of vortex glasses || 
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with anisotropic disorder correlations, which are able to merge the vortex-glass transition 
(equivalent to superroughening) and the Bose-glass transition for j(x,y) = gS(Ax). The 
latter was found at n = 3, as is immediately reproduced by our functional flow equation (^j). 
In fact, this equation holds also in the anisotropic long as the anisotropy of k can 

be neglected. 

The author gratefully acknowledges helpful discussions with J. Krug, T. Nattermann, 
and L.-H. Tang. 
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FIG. 2. Phase-diagram: The free phase, where pinning is irrelevant is bound by 
T/(4irK) n = max(2 — a, 1). For a range of disorder correlations a > a discontinuous 
lock-in transition occurs, whereas for a > 1 a continuous superroughening transition takes place. 
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